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The "true muonium" (/x^/x") and "true tauonium" (r^r~) bound states are not only the heaviest, 
but also the most compact pure QED systems. The rapid weak decay of the r makes the observation 
of true tauonium difficult. However, as we show, the production and study of true muonium is 
possible at modern electron-positron colliders. 
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The possibility of a bound state, denoted here 

as was surely realized not long after the clarifi- 

cation [l[ of the leptonic nature of the muon, since the 
first positronium calculations |2[ and its observation 
occurred in the same era. The term "muonium" for the 
/i+e~ bound state and its first theoretical discussion ap- 
peared in Ref. [J], and the state was discovered soon 
thereafter [s^. However, the first detailed studies 0, 
of (alternately dubbed "true muonium" and 

"dimuonium" [8, 9]) only began as experimental advances 
made its production tenable. Positronium, muonium, 7r/i 
atoms [10|, and more recently even dipositronium [the 
(e+e~)(e+e~) molecule] pjj have been produced and 
studied, but true muonium has not yet been produced. 

The true muonium (/i+/i~) and true tauonium (r+r~) 
[and the much more difficult to produce "mu-tauonium" 
(/i^r^)] bound states are not only the heaviest, but also 
the most compact pure QED systems [the (/i+/i~) Bohr 
radius is 512 fm]. The relatively rapid weak decay of 
the r unfortunately makes the observation and study of 
true tauonium more difficult, as quantified below. In the 
case of true muonium the proposed production mecha- 
nisms include 7r~p ^ (/i+/i~)n [6], 7Z ^ Q, 
eZ^e(/i+/i-)Z [12], Z1Z2 ^ ZiZ2(/i+/i-) [ll (where 
Z indicates a heavy nucleus), direct ji'^ ji~ collisions 0, 



ri^{li'^li )7[M],ande+e -^{fi'^fi ) 



In addition. 



the properties of true muonium have been studied in a 



TABLE I: True fermionium decay times and their ratios. 
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FIG. 1: True muonium level diagram (spacings not to scale). 
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The e+e~ (/i+//~) production mechanism is partic- 
ularly interesting because it contains no hadrons, whose 
concomitant decays would need to be disentangled in the 
reconstruction process. If the beam energies of the col- 
lider are set near threshold ^/s^2m^^ the typical beam 
spread is so large compared to bound-state energy level 
spacings that every nS state is produced, with relative 
probability ~ 1/n^ [i.e., scaling with the (/i+/i~) squared 
wave functions |'0nOo(O)P at the interaction point, r = 
0] and carrying the Bohr binding energy — m^a^/4n^. 
The high-n states are so densely spaced that the to- 
tal cross section is indistinguishable [18] from the rate 
just above threshold, after including the Sommerfeld- 
Schwinger-Sakharov (SSS) threshold enhancement fac- 
tor [19] from Coulomb rescattering. As discussed below 
[Eq. ([2j) and following], the SSS correction na/P can- 
cels the factor of /3, the velocity of either of /i^ in their 
common center-of- momentum (cm.) frame, that arises 
from phase space. 

The spectrum and decay channels for true muonium 
are summarized in Fig.[Tl using well-known quantum me- 
chanical expressions [20] collected in Table H In most 
cases, the spectrum and decay widths of (/i+/i~) mimic 
the spectrum of positronium scaled by the mass ratio 



2 




FIG. 2: The "Fool's ISR" configuration for e+e" ^ (m+m") 
for symmetric beam energies. Tfie angle between the either 
of the and dotted line (z axis) is defined as 0. 



^/i/^e- However, while positronium of course has no 
e+e~ decay channels, (/i+/i~)[n^6'i] ^7* ^ e+e~ is al- 
lowed and has a rate and precision spectroscopy sensitive 
to vacuum polarization corrections via the timelike run- 
ning coupling a{q'^>0). 

Unlike the case of positronium, the con- 
stituents themselves are unstable. However, the fi has an 
exceptionally long lifetime by particle physics standards 
(2.2 /is), meaning that annihilates long before 

its constituents weakly decay, and thus true muonium is 
unique as the heaviest metastable laboratory possible for 
precision QED tests: has a lifetime of 0.602 ps 

in the -^5'o state (decaying to 77) 0,0] and 1.81 ps in the 
^Si state (decaying to e+e~) [6[. 

In principle, the creation of true tauonium (r+r~) is 
also possible; the corresponding and lifetimes are 
35.8 fs and 107 fs, respectively, to be compared with the 
free r lifetime 291 fs (or half this for a system of two 
r's). One sees that the (r+r~) annihilation decay and 
the weak decay of the constituent r's actually compete, 
making (r+r~) not a pure QED system like (e+e~). 

Electron-positron colliders have reached exceptional 
luminosity values, leading to the possibility of detect- 
ing processes with very small branching fractions. The 
original proposal by Moffat [15,] suggested searching for 
X-rays from (/i+/i~) Bohr transitions such as 2P^1S' at 
directions normal to the beam. However, the nS states 
typically decay via annihilation to e+e~ and 77 before 
they can populate longer-lived states. Furthermore, the 
production and rapid decay of a single neutral system 
at rest or moving in the beam line would be difficult to 
detect relative to the continuum QED backgrounds, due 
to a preponderance of noninteracting beam particles and 
synchrotron radiation. 

In this letter we propose two distinct methods for pro- 
ducing a moving true muonium atom in e+e~ collisions. 
In both methods the motion of the atom allows one to 
observe a gap between the production point at the beam 
crossing and its decay to e+e~ or 77 final states. Fur- 
thermore, each given lifetime is enhanced by a relativistic 
dilation factor 7 appropriate to the process. 

In the first method, we utilize an e+e~ collider in which 



at 5 :^ 4m^ carries momentum p = Pg+ -\- Pe- 7^ 0. 
The production point of the (/i+/i~) and its decay point 
are thus spatially displaced along the beam direction. 
Asymmetric e+e~ colliders PEP-II and KEKB have been 
utilized for the BaBar and Belle experiments. How- 
ever, we propose configuring an e+e~ collider to use 
the "Fool's Intersection Storage Ring" (FISR) discussed 
by Bjorken [21] (Fig. [2]) in which the beams are ar- 
ranged to merge at a small angle 26 (bisected by i), so 
that s = (Pe+ + Pe-)^ - 2E+E_(1 - cos2l9) 4m^ and 
the atom moves with momentum pz = (£^+ + E-) cosO. 
For example, for ^ = 5° and equal-energy beams 
E± = 1.212 GeV, the atom has lab- frame momentum 
Pz = 2.415 GeV and 7 = £^iab/2m^ = 11.5. One can 
thus utilize symmetric or asymmetric beams in the GeV 
range colliding at small angles to obtain the cm. energy 
:^ 2m^ for the production of true muonium. 

The gap between the formation of the atom and its 
decay as it propagates should be clearly detectable since 
its path lies in neither beam pipe. The 3^Si state de- 
cays with a 50 ps lifetime, so it moves 1.5 cm of proper 
distance before decaying to e+e~, a length enhanced in 
the lab frame by the 7 factor (to 16.8 cm in the ^ = 5° 
example). One thus can observe the appearance of e+e~ 
events with a ^-dependent set of lifetimes. 

The cross section for continuum muon pair production 
e+e~ ji^ ii~ just above threshold is the Born cross sec- 
tion enhanced by the Sommerfeld-Schwinger-Sakharov 
(SSS) threshold Coulomb resummation factor [1^] S{l3)'. 
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where 



Si/3) = 



Xif3) 



l-exp[-X(/3)] 



(1) 



(2) 



the atomic system produced in e+e 



7* ^ (m+m") 



Here /3 = yl — Am'^/s is the velocity of either of the 

fi^ in their cm. frame, and X{p) =7ray/l — P'^/p. Thus 
the factor of /3 due to phase space is cancelled by the 
SSS factor, so that continuum production occurs even at 
threshold where /3 = 0. For values of \/3\ of order a (as 
in Bohr bound states), we see that the SSS factor effec- 
tively replaces P with ira. Below threshold the entire set 
of ortho-true muonium n^Si C = —1 Bohr bound-state 
resonances with n = 1,2, • is produced, with weights 
~ 1 /n^ and spaced with increasing density according to 
the Bohr energies (v^)n — 2m^— a^m^/4n^. By duality, 
the rates smeared over energies above and below thresh- 
old should be indistinguishable [18]. Thus the total pro- 
duction of bound states in e~^e~ (/i+/i~) relative to 
the e+e~ relativistic lepton pair rate is of order 

R^^na 0:^0. 03. However, in practice the production rate 
is also reduced by the Bohr energy divided by the finite 
width of the beam energies, since only collisions in the 
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FIG. 3: Dominant Feynman diagrams for e^e ). 



energy window 5E are effective in producing 

bound states. For example, if tfie beam resolution is of 
order AE± = 0.01 ~ 1 MeV, the effective R is reduced 
by a further dE/AE± 10~^, leading to a production of 
at ~5xl0~^ of the standard e+e~ ^ £^£~ rate. 

Since the state in the FISR method is pro- 

duced through a single C = — photon, n^Si states 
(and not n^So) are produced, which decay almost al- 
ways to pairs, as illustrated in Fig. [2l Note that this 
holds even for radiative transitions through the sequences 
n"^Si n' P n^Si^ since the intermediate P states 
do not annihilate. 

The bound states, once produced, can in prin- 

ciple be studied by exposure to 0(ps) laser or microwave 
bursts, or dissociated into free /i^ by passing through a 
foil. Because of the novel kinematics of the FISR, the 
true muonium state can be produced within a laser cav- 
ity. For example, an intense laser source can conceivably 
excite the nS state of the atom to a P state before the 
former's annihilation decay. A 2P state produced in this 
way has a lifetime of (15.4 ps)x7. In principle this al- 
lows precision spectroscopy of true muonium, including 
measurements of the 2P-2S and other splittings. Laser 
spectroscopy of {jaZ) atoms is reviewed in Ref. [22] . 

In this letter we also propose a second production 
mechanism, e+e~ 7(/i+/i~), which can be used for 
high-energy colliding beams with conventional configu- 
ration. It has the advantage that the production rate is 
independent of beam resolution, and removes the (/i+/i~) 
completely from the beam line since the atom recoils 
against a co-produced hard 7. While the production of 
the real 7 costs an additional factor of a in the rate, the 
kinematics is exceptionally clean: Since the process is 
quasi-two-body, the 7 is nearly monochromatic [neglect- 
ing the (/i+/i~) binding energy] as a function of the total 
cm. beam energy ^/s^ E^ = {s — 4m^)/2v^. Further- 



more, the {ji'^ lifetime is enhanced by the dilation 
factor 7 = (5 + Am'jJ/ 4:171^^. The dominant Feynman 
diagrams for e+e~ -^'y{i~^i~) are shown in Fig. [3l If E^ 



is large compared to m^, and its angle from the beam is 
large, one also can have significant radiation from the ja 
lines. In fact, if the hard 7 is emitted by a /i, the true 
muonium state is formed in the para n^So C = state, 
which leads to two-photon annihilation decays. In this 
case both e+e~ and 77 final states appear, accompanied 
by a decay gap. 

The 0{a^) Born amplitude for the process e+e~ 
^ji^ li~ (free /i's) in the kinematic regime ml/ m^/s<Cl 
was first computed long ago in Ref. [23]. More recently, 
a related collaboration [24] specialized this calculation 
to precisely the desired kinematics: The invariant mass 
square si of the ji^ pair is assumed small compared to 
the total cm. squared energy s. In this case, the Born 
differential cross section for the diagrams in Fig. [3] is 



d(j 



SSi{l — c^) 



(25 + 1 — 2x_x+) dx- dcdsi 



(3) 



where c is defined as the cosine of the angle between 
the e~ and /i~ [and hence also the (/i+/i~) atom], 6 = 
w?^/ si :^ I for {ji^ bound states, x± are the fractions 
of the half of the cm. beam energy E±/{^/2) that is 
carried by /i^ (the other half being carried by the 7) so 
that x_ = 1, and the range of x± is given by 



l{l-(3')<x±<\{l 



■/?'), 



(4) 



where = \^1 — 4S is the velocity of either of the /i^ in 
their cm. frame. In addition, the hard photon momen- 
tum makes an angle 0^ with the initial e~ that is assumed 
to lie outside of narrow cones of opening angle Oq around 
the beam axis, ^0 < ^7 < tt — 6>o, where 2m^/v^<C^o 1- 
Note that the 7 and (/i+/i~) are back-to-back, 0^ = 7t—0. 

The differential cross section in Eq. (|3]) becomes singu- 
lar when the 7 [and hence also the (/i+/i~)] is collinear 
with the beam. For the purpose of our cross section esti- 
mates, we integrate c over the range excluding the beam 
cone, cG [— co,co], where Co = cos6>o- Using also Eq. (|4|) 
to integrate over one obtains 



da 
dsi 



213' 



In 



I+Cq 
1 - Co 



Co 



a 

SSi 



(5) 



The factor jS' ^ indicating that the cross section vanishes 
at (1^ threshold, arises simply from 3-body phase space. 

Equations (|3]) and ([5]) describe a process in which the 
fi^ pair carry an invariant mass si small compared to 
5, but are not necessary bound together. In order to 
compute the cross section for such a process, one must 
again include the SSS threshold Coulomb resummation 
factor [19]. Here the (5' in Eq. ([5j) refers to the (con- 
tinuous) velocity of each of a free /i^ pair in its cm. 
frame, whereas (3' in the bound-state formalism refers 
to the (quantized) velocity of each particle within their 
Coulomb potential well. Nevertheless, as argued in the 
previous case, by duality the same cross section formulae 
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still hold in the bound-state regime if the SSS factor is 
taken into account and the weights of the discrete transi- 
tions are properly included and smeared over the allowed 
energy range for bound states [18]. One then obtains 
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The relevant range of dsi is just that where bound Bohr 
states occur, which begin at energy a^m^/A below the 
pair creation threshold si= 4m^, and thus give rise to 
dsi'^m^^a^. Thus one obtains 
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(7) 



The angular factor is again singular for co = ±l, varying 
from zero at ^o = 7r/2, to unity near 7r/4, to over 7 at 2°. 

Note that /3 = ^l — Am?^/s differs from (5' used for 

e+e~ for processes e+e~ 7(/i+/i~) with 

the same value of s for which Eq. (|3]) and following are 
applicable, recall that <C s and hence 1. The ratio 
cr[e+e~ ^7(/i+/i~)]/cr(e+e~ at the same value 

of s is therefore just a number close to unity (e.^., 2.66 for 
Oo = 2 degrees) times a^. While this O(10~^) suppres- 
sion may seem overwhelming, it is within the capabili- 
ties of modern e+e~ facilities. For example, the BEPCII 
peak luminosity will be 10^^cm~^s~^ at a cm. energy of 
3.78 GeV, but varying between 2 and 4.6 GeV [25]. At 
2 GeV about 5 events e+e~ ^7(/i+/i~) occur per year 
of run time, and the yield increases with 1/s. On the 
other hand, for smaller values of s the dilation factor 7 
for becomes shorter, thus diminishing its lifetime 

and hence track length. 

The production is much more prominent if one per- 
forms a cut on si values near the /i^ threshold 4m^. In 
that case one should compare Eq. ([5]) to the derivative 
da{e~^e~ fi~^fi~)/ds [which is not actually a differen- 
tial cross section but rather the difference of cr(e+e~ 

between bins at cm. squared energies s and s-\-ds]. 
Then the relative suppression is only O(a^), one a arising 
from the extra photon and one from the SSS factor. 

Between the two proposals presented here, e+e~ 

with beams merging at a small crossing angle, 
and the rarer process e+e~ 7(/i+/i~) that can access 
both ortho and para states with conventional beam kine- 
matics, the discovery and observation of the true muo- 
nium atom appears to be well within current 

experimental capabilities. 
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